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Abstract 

We investigate two methods of obtaining exactly solvable potentials with 
analytic forms. 



PACS numbers:03.65.Ge,11.30.Pb 



There are two methods of obtaining exactly solvable potentials in quantum me- 
chanics. The first method was developed by applying the technique of supersymmetry 
(SUSY) to the Schrodinger equation and obtain two potentials with almost identical 
spectra. The two potentials can be considered to be superpartners of each other. 
It has been shown |jl|, |] that if the two partners happen to be related by a simple 
relationship called shape invariance, the energy eigenvalues of the potential can be 
solved exactly. All the known solvable potentials with closed analytic forms can be 
shown to be shape invariant. In the literature, there are other solvable potentials 
that have not been shown to be shape invariant, however, they exist only complex 
numerical forms which we shall not consider in this article. 

A second interesting method of obtaining solvable potential was proposed by 
Klein and Li^l based on some special quantum commutation relationships. Li|^ has 
recently worked out the most general potential that can be obtained this way. 

In this paper, we first investigate the relationship between the two approaches. 
We will show that the general solutions obtained by Li are special cases of the general 
solutions that can be obtained by solving shape invariance condition. Therefore, the 
solving shape invariance condition remains the most general method of obtaining 
exactly solvable potentials which can be given in analytic form. Unfortunately, there 
is no general method for getting analytic solutions of shape invariance condition. 
The best one can do seems to be starting from a guessing ansatz with an unknown 
function. The shape invariance is then enforced by demanding that the function 
satisfies an ordinary differential equation. The ansatz allows one to turn the difficult 
shape invariance condition into a problem of solving differential equation. It also 
allows one to associate each ansatz as defining a particular class of solutions. While 
the solutions obtained by Li correspond to those defined by a particular ansatz, 
in the literature, the largest classes of analytic solutions of the shape invariance 
condition was provided by Gendenshtein He proposed three ansatzs which define 
three classes of solutions which seem to cover all the known analytic exactly solvable 
potentials. We, therefore, proceed to solve the differential equations corresponding 
to these ansatzs and the energy eigenvalues of these three classes of solutions of 
shape invariance condition. In the process, we also demonstrate that all the known 
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solutions in the hterature are special cases of the one provided by Gendenshtein. 

The three ansatzs of Gendenshtein and the forth one corresponding to Li's so- 
lutions also demonstrated that each potential can be simultaneously represented in 
many classes. It is also clear that one can in principle continue to invent more ansatz 
without obtaining any new potentials. Therefore, one is immediately faced with the 
intriguing question of what is the better way to classify these solvable potentials 
than using the ansatzs. One method seems to be starting from the n dependence 
of its spectrum. It is noted that only a few simple general forms of n dependence 
are allowed for the known exactly solvable potentials. We shall make some com- 
ments in this direction later and mention some interesting unsolved problems in the 
conclusion. 

We first briefly review the procedure of constructing creation and annihilation op- 
erators of simple harmonic oscillator, now we can apply this method to the Schrodinger 
equation with an arbitrary potential V{x). 

Assuming that we have already shifted the potential by a constant so that the 
ground state energy becomes zero, we shall denote this Hamiltonian and its potential 
with (— ) subscript and denote ip'^'' as its eigenf unctions. That is, Vl(a;) = V{x), 
4'\x)=Mx), = 0, with 

H-4-^ = {-l-^,+V4x))4-^ = 0. (2) 

The above equation, Eq.(|^) is identical to 



1 r/2 7//-)" 
--^(-^^-V (3) 



Then, the general creation and annihilation operators are: 



V2^ dx 4")^' 



1 ^ d 



The Hamiltonian can be written as H_ = A. Now we define a new Hamiltonian 
if+ = AA^ . which can be written as if+ = ~\-^ ^ ^+{x). The corresponding 



3 



potential in the new Hamiltonian is: 

V,^-V_^iP^^V_-^^nii-> (5) 

The V^, are called supersymmetric partner potentials. 

One can define a new function, W{x), called superpotential. 

Solving the differential equation we get 

4')(a;) = exp{-V2 J'' W{x')dx'), (7) 

and = —^-^ + W{x), A = -^-^ + W{x). The two partner potentials are made 
from superpotential W{x), 

A'^ and A do not commute with each other and satisfy [A, A"*"] = \/2W{x)' . To 
compare the spectra of these two partners, denote as the eigenfunction of if_ 
and ■j/'*^^-' as that of Then, Ail)^~^ is an eigenfunction of iJ+ because 

Similarly, A^^/'^"'"^ is an eig enfunction of H^, H^{A+i)^+^) = There- 
fore A, A'^ can transform an eigenfunction of one potential into a eigenfunction of 
it's partner with the same energy. However, note that the ground state of V- is 
annihilated by A and have no partner state of V+. Therefore the ground state of V- 
has no superpartner. 

Relations of energy eigenvalues and wavefunctions are: 

^1+^ = ^i=M;\, n = 0,1,2,..., 

^it = -7^A^i^i^\ n = 0,1,2,.... (10) 

In 1983, Gendenshtein [||] suggested that all exactly solvable potentials are "shape 
invariant" ( he used "form invariant" instead ). Shape invariance means that the 
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superpartner potentials have the same x— dependence modulo the changes in val- 
ues of a set of parameters that define the potential. Mathematically, it means the 
superpartner of V-{x; ao), V+{x; oq) can be written as 

V+{x; ao) = V_{x; ai) + R{ai) (11) 

and Oi = /(oq), where Oq are the set of parameters in V^, and / is the transformation 
function that maps into Oi. The remainder R{ai) is independent of x. 

To show that the eigenvalues can be obtained easily from the above condition, we 
construct a series of Hamiltonian H'-''\ A; = 0, 1, 2, . . ., with H^^^ = H^, H^^^ = H+, 

H^'^ = + Vl(x; a,) + ^ i?(a,), (12) 

where = /'^(ao), i.e., / mapped s times. Furthermore, we have From Eq.(0), we 
see that the ground state of H^'^^ has the energy: 

k 



Eit^ = Y.Rias). (13) 

s=l 



Since (n + l)th energy eigenvalue of H^^^ (=i/_), whose ground state energy is 
zero, is coincident with the ground state energy of Hamiltonian H^'^\ the complete 
eigenvalues of H are: 

Ei^ = X: Ria,) , Et^ = 0. (14) 

k=l 

If a potential is shape invariant, we can also get the bound state wavefunctions 
■?/^^~^ easily. This is because A and A'^ can link up the wavefunctions of the super- 
partners with the same energy. Starting from H^~\ it's ground state iljQ~\x;an) 
corresponds to the first excited state tpl ^(x;a„_i) of H^Ji- In the same manner, 
eventually, it will correspond to the nth state of i/g Recalling Eq.(|l^), one ob- 
tains 

ipl^^ oc A+(x; ao)A^{x; ai) . . . A+(x; an-i)i'o~\x; an). (15) 

To compare the supersymmetry approach with the second approach proposed by 
Klein and Li^, |^, we shall first give a brief review. Consider an one dimensional 
quantum system with the Hamiltonian 
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For an arbitrary function of position, / = /(x), one can easily derive 

[[/, HIH] = -if'H + Hf") + 2f"V + f'V - (17) 

To obtain solvable potential Klein and Li proposed to impose a so-called linear double 
commutator relation demanding that for some / and V, the right hand side become 
a linear functional of /. In order to make the operator equation linear in /, one 
imposes 

2/'V + /V' = «/ + /?, (18) 

and 

f" = + (19) 
Eq. ([l8|) relates to / and can be solved to give the potential V{x) 

Vi^^-'^^^- (20) 

Eq. ([l9D can be solved for / as 

f = ax^ + hx + c (/i = 0), (21) 

or 

/ = AeV^^ + 5e-v^^ + C (yU^O). (22) 

For /i = 0, the resulting potentials are just those for one- or three-dimensional 
harmonic oscillator problems. For yU 7^ 0, one can get the general potential 

which are special cases of Morse and Poschl- Teller potentials. Li also solved the 
energy eigenvalues by the method of Heisenberg's matrix mechanics 0] under some 
general assumptions. The eigenvalue is proportional to n^, where n is a quantum 
number. Also, from Eqs. ([T9|) and (0), the corresponding superpotential, with 
linear double commutator relation, can be written as[^] 

Here we shall first show that the superpotential is indeed shape invariance and then 
solve the eigenvalues using the supersymmetric methods outlined earlier. 



6 



Given, Eqs.(24) and ([T9|) , one can easily check the shape invariance condition by 



working out the potentials V± as 

T ^,)f + {ajb + bja T^))f + bjb T ^) a 

- Ur ^ 7=2- ^''^ 

Therefore, V± are related by shape invariance condition 

V+{a,b-x)=V4a--^,b-^;x) + V2a-i^. (26) 
The energy eigenvalues can be straight-forwardly worked out to be 

2 

En = V2na - (27) 

where the ground state energy has been shifted to zero. SUSY method and Li's 
method gives the same results as expected. 

In the SUSY approach, shape invariance requires that the superpotential W{x) 
satisfies the functional differential equation: 

W\a, x) + -^W'{a, x) = W\ai, x) - -^W'{ai, x) + R{ai), (28) 
V 2 v 2 

where Cq represents a set of parameters, called the shape invariance relation(SIR). 

The values of the set, Oi, on the right-hand side, depend on the value of Oq, that is, 

ai = /(ao) for some function /. The only known way of solving SIR is to impose an 

educated ansatz to turn it into a differential equation. In this direction, Gendenshtein 

|jl|] had proposed three ansatzs which provide some general classes of solutions. There 

is no claim that the three classes should encompass the most general solutions of SIR. 

However, it is interesting to note that all the known solutions can be shown to be 

special cases of one of the three classes as we shall demonstrate later. Also, the three 

classes are not mutually exclusive. That is, some solutions can be represented in 

more than one classes. 

We shall discuss the three classes in order, solve the corresponding nonlinear 

differential equations and work out the potentials and energy eigenvalues for each 

class. The three classes can be described as follow: 

(I). Class one: 

The ansatz for the superpotential is of the form 

W = afi + b. (29) 
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The SIR then requires /i to satisfy 



f[=pfi+qfi+r. (30) 

V± can be worked out to be 

V± = a{a ± ^)f^ + 2a{b ± ^) A + b'±^, (31) 
with the parameters transforming as 

np 



9 I "^Pg 



a6 ^°72 + ^ 



R = bl-bl + —^{ao + ai). 
v2 

For the case when p 0, the eigenvalues can be worked out to be 

n 
k=l 

k=i v2 
Ordering the terms by power of n, En can be written as 

. , _ _ , _ , _ , ,| ,34) 

For the case when p — 0, the transformation of parameter and the eigenvalues can 

be obtained as 

bn = b+^, (35) 



R = bl-bl + V2ar, 



En = ^ -Rfe = V2{ar - bq)n - —n^. 
k=i ^ 

The known exactly solvable potentials in this class and their corresponding param- 
eters are hsted in Table 1. 
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Table 1: Class One W = af{x) + b ; f = pf + qf + r 



Potential 


W{x) 


/(^) 


a 


b 


P 


Q 


r 


Shifted Oscillator 


^UJX-b 


X 




-b 








1 


Coulomb 


[T e2 l+l 
V 2 l+l V2r 




V2 


[T e2 
V 2 Z+1 


-1 








Morse 


A - Se-"=" 




-B 


A 





—a 





Rosen-Morse 


A tann ax + 


tanh ax 


A 


D 


—a 





a 


Eckart 


—A coth ar + ^ 
{B > A") 


coth ar 


-A 


B 
A 


—a 





a 



The most general solutions for function / in this case can be summarized as 

f^rx + c ,(p = g = 0); (36) 

/ = ^— ,(p7^0,? = r = 0); (37) 

px + c 

/ = ce^^-^ ,(p = 0,5 7^0); (38) 

V4pr - V4pr - q 
/= ^- tan( ^ + ^)~^ ,(P7^ 0,g or r 7^ 0), (39) 

where c is an integration constant. 
(11). Class two: 

The superpotential in this ansatz is assumed to be of the form 



J2 



(40) 



V± can be worked out to be 

V2 



V+ = a{a ± -^)f^ + b{b T ^)/2-' + 2a6 ± - Pb), (41) 
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with the parameters transforming as 

np 

Rn = 2V^(ag - hp) + 2(2n - l)pq. (42) 

The eigenvalues are 

En = 2V2{aq - hp)n + 2pgnl (43) 

The well-known examples in the class and their corresponding parameters are listed 
in Table 2 for illustration. 



Table 2: Class Two W ^ af + ^ ; f ^ pf + q 



Potential 


ir(.r) 


/(•O 


(I 




P 


(1 


3-D Oscillator 




r 




i+i 





1 


Poschl-Teller I 


A tan ax — B cot ax 


tan ax 


>1 


-B 


a 


a 




{0 <ax < I) 












Poschl-Teller II 


A tanh ar — B coth ar 


tanh ar 




-B 


—a 


a 




(5<A) 













The general solutions for function / in this case can be written as 

(44) 

(p^ 0,^7^0), (45) 

= 0), (46) 

where c is an integration constant. 
(III). Class three: 

The superpotential is assumed to be of the form 

^ Q + bsjpfi + q 
h 

fs = \/pfi + q- (47) 



f^qx + c (p = 0), 

jq 

f = J - ta.n( ^/pqx + c) 
\ P 



f 



1 

px + c 



{p^O,q 
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V± can be worked out to be 



1 . o , „ 1. „ 1 



V± = jiia' + bq{b T y 2) + 2«(^ T + + b\ 



with the parameters transforming as 



The eigenvalues are 



2 

= V2npb - ^ 



The well-known exactly solvable potentials in this class are listed in Table 3. 
The general solutions for function / are 

/ = ^x + c ,(p = 0); 
f^ce^^ ,(p^O,g = 0); 

c is an integration constant. 



Table 3: Class Three W^{a + by/pP + q) // = y/pP + q 



Potential 


ir(,r) 


,f(.r) 


(1 


h 


P 


(1 


Morse 


A - £e-"^ 




-B 


A 
a 


a^ 







A tanh ctx -|- S sechctx 


cosh ax 


B 


A 
a 


a^ 


-a2 


Rosen- Morse 


A coth ar — B eschar 

{A<B) 


sinh ar 


-B 


A 
a 


a^ 


a^ 


Eckart 


—A cot ax -|- S CSC ax 
(0 < ax < TT, A > S) 


sin ax 


B 


_A 

a 


-a2 


a2 
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The tables showed that all known exact solvable potentials with analytic forms 
can be put into one of these three classes. Whatever has spectrum \jr? can be 
classified as the first case. However, only the class two solutions can produce three 
dimensional oscillator, type-I Poshl- Teller or type-II Poshl- Teller potentials. But 
two potentials, Morse and Eckart, can be considered both as case one and as case 
three. As mentioned before, Li's results are just the special case of shape invariance 
solutions. However it cannot be so easily fit into one of the three classes given by 
Gendenshtein. In fact one can consider Eq.(|T9D and Eg . (plD to be the equations that 
define a fourth ansatz for solutions of SIR. The set of solutions overlaps with those of 
the other three classes provided by Gendenshtein but does not generate a new one. 

From this point of view, it is clear that the ansatz does not provide a very precise 
classification of the solutions of shape invariance condition. It is not too hard to 
propose new ansatz, however, it is much harder to generate new solutions. Typically 
one can not be sure whether an ansatz generates any new solution or not until they 
are solved completely. Therefore it seems that a better classifying solutions may be 
to use the energy spectrum and its quantum number dependence instead. 

In conclusion, we have discussed the two methods of obtaining exactly solvable 
potentials in quantum mechanics. One of them requires the shape invariance be- 
tween the superpartners of the potentials in the supersymmetric formulation. The 
other one imposes a so-called "linear double commutator relations" . We have shown 
that the second case only produces solutions which are special solutions of the first 
approach. From this point of view, the shape invariance approach is still the most 
general method of producing the analytic, exactly solvable potentials. We also ar- 
gued that the n dependence of the energy spectrum may be a better way of telling 
the difference between different classes of potentials. In addition, we work out the 
energy eigenvalues of the most general classes of potentials in the literature. 

Unfortunately, it is still not possible to obtain the most general solutions to SIR. 
From the table of 0, one observes that all the known shape invariant potentials 
have basically only one parameter changing under the SUSY transformation. In 
particular, for Poschl- Teller I potentials, only the combination A + i? is changing. 
The parameter A — B vs, invariant. For Poschl- Teller potentials II, on the other 
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hand, the combination A — B is changing and the parameter A + B is invariant. It is 
not surprising because as long as the transformation property of the parameters are 
hnear for a proper choice of parameters, one can always make linear combinations 
such that only one of them is changing during the transformation. Since, in general, 
there is no reason for the transformation of the parameters to be linear, one would 
expect to have a lot more interesting solutions of shape invariance condition waiting 
to be discovered. 

Another interesting observation is that all the shape invariant solutions has the 
spectrum of one of the following the forms (modulo a constant that sets the ground 
state energy to zero): (1) an, (a > 0); (2) — (^^ip^, (b > 0); (3) ±(a + bny or their 
linear combination. The harmonic oscillator is an example for the first form. The 
hydrogen atom is an example for the second form and, the square well is the simplest 
example for the third form. Also note that the general spectra of all three classes of 
solutions suggested by Gendenshtein are all linear combinations of the above three 
forms. One may wonder if there is something fundamental about these kinds of 
spectra that made them represent the spectra of all the shape invariant, exactly 
solvable, potentials. 

Finally, regarding the approach proposed by Klein and Li, since Li's solutions in 
0] produce only part of the solutions of SIR, it suggests that it may be possible to 
generalize their approach within its' framework. We have made some attempts in 
this direction, however, so far without success. 

This work was supported in part by National Science Council of Republic of 
China grants No. NSC 84-2112-M-007-042 and No. NSC 84-21 12-M-007-016. 
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